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ABSTRACT 



We present a kinematic analysis of the globular cluster (GC) system in the 
giant elliptical galaxy (gE) M60 in the Virgo cluster. Using the photometric 
and spectroscopic database of 121 GCs (83 blue GCs and 38 red GCs), we have 
investigated the kinematics of the GC system. We have found that the M60 GC 
system shows a significant overall rotation. The rotation amplitude of the blue 
GCs is slightly smaller than or similar to that of the red GCs, and their angles of 
rotation axes are similar. The velocity dispersions about the mean velocity and 
about the best fit rotation curve for the red GCs are marginally larger than those 
for the blue GCs. Comparison of observed stellar and GC velocity dispersion 
profiles with those calculated from the stellar mass profile shows that the mass-to- 
light ratio should be increased as the galactocentric distance increases, indicating 
the existence of an extended dark matter halo. The entire sample of GCs in M60 
is found to have a tangentially biased velocity ellipsoid unlike the GC systems 
in other gEs. Two subsamples appear to have different velocity ellipsoids. The 
blue GC system has a modest tangentially biased velocity ellipsoid, while the 
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red GC system has a modest radially biased or an isotropic velocity ellipsoid. 
From the comparison of the kinematic properties of the M60 GC system to those 
of other gEs (M87, M49, NGC 1399, NGC 5128, and NGC 4636), it is found 
that the velocity dispersion of the blue GC system is similar to or larger than 
that of the red GC system except for M60, and the rotation of the GC system is 
not negligible. The entire sample of each GC system shows an isotropic velocity 
ellipsoid except for M60, while the subsamples show diverse velocity ellipsoids. 
We discuss the implication of these results for the formation models of the GC 
system in gEs. 

Subject headings: galaxies: clusters: general — galaxies: individual (M60) — 
galaxies: kinematics and dynamics — galaxies: star clusters 



1. Introduction 



Globular Clusters (GCs) have been long recognized as an important tracer to under- 
stand the formation and evolution of galaxies. To solve the mystery of galaxy formation and 
evolution, several photometric properties of th e GCs such as color distribution, spatial struc- 
ture, and luminosity function have been used (lLedl2003l : IWest et al.ll2004l : iBrodie fc Strader 
20061 ). However, it is still difficult to test the predictions of the formation models of galaxies 
and their GC systems using photometric data alone. 

Recently, with an aid of large telescopes with apertures larger than 4 meter, a large 



statistically mean ingful kinematic studies (e.g. 



Cote et al. 


2001. 


2003: 


Richtler et al. 


2004 



Peng et al.l 120041 ). Due to their wide spatial distribution, brightness, and compact size, 
GCs are useful test particles to trace the gravitational potential of their host galaxies, and 
especially of the dark matter halo beyond several effective radii of the galaxy. Therefore, 
a kinematic study of the GC system enables us to estimate the global mass distribution 
of their host galaxy, or to constrain the orbital properties of GCs using an independently 
determined mass profile (e.g., from X-ray emission) of the galaxy. Moreover, the kinematic 
difference between GC subpopulations (blue and red GCs) can be used as an observational 
constraint on the galaxy formation model. 

To date there are six giant elliptical galaxie s (gEs) for which the kinematics of the GC 
system have been stud ied to o ur knowledge: M87 (ICohen fc Rvzhovlll997l : iKissler-Patig fc Gebhardt 



1998 



1998 



Minniti et al. 



Cote et al.ll200lh. M49 JZeof et al. 



2000 



Cote et al.ll2003r). NGC 1399 (iKissler-Patig et al. 



19981 : iKissler-Patig et al.lll999l : lRichtler et al.ll2004l ). NGC 5128 flPeng et al. 
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2004 : IWoodlev et al.l 120071 ) . NGC 4636 flSchuberth et all l2006h . and M60 flBridges et a 
20061 ). M87, cD galaxy of the Virgo cluster, was studied recently by ICote et al.l (l200ll ) 
using ~ 280 velocity data of GCs. The velocity dispersions of the blue and red GCs were 
found to be ~410 km s^^ and ~390 km s~^, respectively. Both the blue and red GCs appear 
to rotate around the photometric minor axis with a similar rotation amplitude of ~ 160 km 
s^^ when averaged over the whole system, while the blue GCs appear to rotate around the 
photometric major axis inside a radius of ~ 16 kpc. The entire GC system has an isotropic 
velocity ellipsoid, while the b l ue an d red GC systems show tangentially and radially biased 
one, respectively. ICote et al.l (120031 ) studied the GC kinematics of M49, the brightest mem- 
ber of the Virgo cluster, using ~ 260 velocity data of GCs. They found a larger velocity 
dispersion of the blue GCs (~350 km s~^) than that of the red GCs (~270 km s~^). The 
blue GCs rotate roughly around the photometric axis of M49 with a rotation amplitude of 
~ 100 — 150 km s~^, while the red GCs show some evidences for weak rotation (~ 50 km s~^) 
around the same axis but a different direction. The entire, blue, and red GC systems appear 
to have an isotropic velocity el lipsoid. NGC 1399, cD galaxy of the Fornax cluster, was 
studied by iRichtler et al.l (120041 ) with the largest number (~ 470) of GC velocity data. The 
velocity dispersions of the blue and red GCs were estimated to be ~291 km s~^ and ~255 
km s~^, respectively. No significant rotations were found for both the blue and red GCs, 
while there is a weak signature of rotation for blue GCs beyond 6'. The velocity anisotropies 
for both th e blue and red GC s are consistent with isotropic orbits. Using ~ 220 GC data of 
NGC 5128, iPeng et al.l (120041 ) found that the red GCs ex hibit a significant rota tion, while 
the blue GCs do not show a clear hint of rotation. Later, IWoodley et al.l (120071 ). using 340 
GC data for NGC 5128, showed that the rotation a mplitude and the veloci ty dispersion for 
the subsamples are quite similar. For NGC 4636, ISchuberth et al.l (120061 ) found that the 
velocity dispersions of the blue and red GCs are not different, but the rotation of the red 
GCs is stronger than that of the blue GCs. In summary, all these GC systems show diverse 
GC kinematics in velocity dispersion, rotation, and their radial variation (see Section H?T]) . 
making it difficult to draw any strong conclusion on the uniform formation history of the 
GC system in gEs. 



M60 (NGC 4649) is a giant elliptical galaxy in the Virgo cluster, slightly less lumi- 
nous than M87 and M49. M60 has a nearby companion Sc galaxy, NGC 4647, located 
at 2. '5 from the center of M60. While there were num erous photometric s t udies of the 
GC system in M60 based on the ground-based images (ICouture et al.l Il99ll : iHarris et al. 



19911 : lAshmari fc Zepi 
scope images ( NeilsenI ' 



Strader et al. 



200 



i 



1998 



1999 



Forbes et al. 



200 



i 



Kundu fc Whitmore 



Mieske et al 



system (jPierce et al 



Lee et al 



2001 



20071) and Hubble Space Tele- 



Larsen et al.ll200ll : IPeng et al.ll2006 



200f 



^ , there w ere few spectroscopic studies of the M60 GC 



20061 ). iPierce et al.l (120061 ) published for the first time the spectroscopic 
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observational results for 38 GCs (16 blue and 22 red GCs) in M60 using Gemini/GMOS. 
They found no obvious signs of a recent starburst, intera ction or merge r by es timating the 
ages and metallicities from the spectra of M60 GCs. Later. [Bridges et al.l (120061) i nvestigated 
the kinematics of M60 GC system using the velocity data of iPierce et al.l (120061 ). They re- 
ported that the velocity dispersion of the blue GCs is smaller than that of the red GCs unlike 
the cases of other gEs. They found no hint of rotation in the GC system of M60. Further- 
more, orbital distribution of the GC system is found to be close to isotropic inside the radius 
of 100", but becomes tangentially biased beyond this radius. However, due to their small 
number [N ~ 38) of GC velocity data with limited radial and azimuthal coverage, it was 
difficult to derive the kinematic properties up to several effective radii, and to distinguish 
the kinematic difference between the blue and red GCs. 

We have carried out a photometric study of t he GC system i n M60 using deep wide- field 
CTi images obtained at the KPNO 4m telescope (ILee et al.ll2007l ). and a spe ctroscopic study 



using the spectra obtained at the Canada-France-Hawaii Telescope (CFHT) (ILee et al.ll2006 



Paper I). In this paper, we present the results of kinematic study of the GC system in M60 
using the velocity data of 121 GCs selected in Paper I. Section 2 gives a brief description 
of the data used in this analysis, and the kinematic properties are derived in Section 3. 
From the comparison of kinematic properties of the M60 GC system to those of other gEs, 
we discuss our results regarding the GC formation models in Section 4. A summ ary of this 



study is given in the final section. We adopt a distance to M60, 17.3 Mpc, given by lMei et al. 



( 120071 ) based on a surface brightness fiuctuation method, and the corresponding scale for one 
arcmin is 5.032 kpc. 



Data 



We used the spectroscopic data of GCs given in Paper I, which describes the details of 
the spectroscopic observation, data reduction, and the data set. Here we only give a brief 
summary of the data set of M60 GCs. 

We selected GC candidates in deep, wide- field Washington C and Ti images (16' x 16') 



obtained at KPNO 4m telescope (ILee et al. 



in th e HST/WFPC2 archive data (iNeilsen 



2007) and F555W (V) and F814W (I) images 



1999: Kundu fc Whitmore 2001: Larsen et al. 



200ll ). Spectroscopic observations were made using the Multi Object Spectrograph (MOS) 
at the 3.6 m CFHT in February 2002 and in May 2003 for 165 GC candidates with 19 < 
Ti < 22 mag and 1.0 < (C — Ti) < 2.4. We determined the radial velocities of GC 
candidates by cross-correlating the candidate spectrum with that of three Galactic GCs. 
Among 165 GC candidates, we could extract the spectra and determine the radial velocities 
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for 111 objects, except for 54 objects due to the poor quality of their spectra. We increased 



the nu mber of GCs by combining ou r data of radial ve locities with those of iPierce et al. 



( 120061 ) ■ The radial velocities of GCs in lPierce et al.l (120061 ) were transformed into our velocity 
system using equation (1) in Paper I, and the transformed velocities were used for further 
anal ysis. Of the entire spectroscopic sample of GC candidates (110 from Paper I and 38 
from lPierce et al.l 120061 ). 121 genuine GCs were selected in Paper I using a criteria of radial 
velocities (500 < Vp < 1600 km s"^) and {C - Ti) colors (1.0 < C - Ti < 2.4 mag). There 
are 83 blue GCs with 1.0 < (C - Ti) < 1.7 and 38 red GCs with 1.7 < (C - Ti) < 2.4 
i n the total sample. Foreground reddening toward M60 is very small, E{B — V) = 0.026 
dSchlegel et aLlll998h . corresponding to E{C -Ti) = l.966E{B-V) = 0.051, A{Ti) = 0.071, 
and A{V) = 0.088. 

In Figured], we show the spatial distribution of 121 GCs in M60 with measured velocities. 
Two large dotted ellipses represen t isophotes at the 25.0 B-m ag arcsec"^ of M60 (larger 
one) and NGC 4647 (smaller one) (jde Vaucouleurs et al.lll99ll ). It is worth noting that the 
majority of high velocity GCs (open symbols) with larger velocities than the systemic velocity 
of M60 (Vfrr, E = 1056 ± 64 km s-\ Paper I), appear to be located to the north-west of M60. 
On the other hand, the majority of low velocity GCs (filled symbols) with smaller velocities 
than the systemic velocity of M60 appear to be located to the south-east of M60. It might 
be related to the overall rotation of the M60 GC system around the minor axis. However, a 
careful analysis is needed for further dynamical investigation due to the non-uniform spatial 
coverage of observed genuine GCs around M60. 



3. Results 

Using the master catalog of 121 GCs in M60 (Paper I), we have investigated the kine- 
matic properties of the M60 GC system: the rotation amplitude, the position angle of the 
rotation axis, the mean line-of-sight velocity, the projected velocity dispersion, and the ve- 
locity eUipsoid. 



3.1. Rotation of the Globular Cluster System 



A detailed description of the relation betwe en the intrinsic rotational velocity field and 

(1200 ll ). In summary, if we assume that the GC 



the projected one is given in iCote et al 



^Vg^i = 1117 ±6 km in lConzalezI |l993h . 
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system is spherically symmetric with an intrinsic angular velocity field stratified on spheres, 
and that the GC rotation axis lies in the plane of the sky, then we can find that radial 
velocities of GCs depend sinusoidally on the azimuthal angles. For the M60 GC system, 
the assumption of spherical symmetry of the GC s ystem is reasona ble due to the modest 
projected ellipticity (effective ellipticity, Ces = 0-21. iLee et al.l 120071 ). Therefore, we fit the 
observed hne-of-sight velocities (vp) of the GCs with the function. 



vp{e) 



Jsjs 



+ (QR) sin(e - Go) 



'1^ 



where B is the projected position angle of GCs relative to the galaxy center measured 
from north to east, 6o is the projected position angle of the rotation axis of the GC system, 
R is the projected galactocentric distance, (QR) is the rotation amplitude, and fsys is the 
systemic velocity of the GC system. 

In Figure[2], we plot the radial velocities of GCs with measured uncertainties as a function 
of position angle for entire 121 GCs (top panel), 83 blue GCs (middle panel), and 38 red 
GCs (bottom panel). The best fit rotation curve of equation ([T]) for each sample is overlaid. 
The fitting was done using an error-weighted, nonlinear fit of equation ([1]) with Vsys as a fixed 
value of M60 recession velocity (fgai = 1056 ± 64 km s~^) rather than a free parameter for 
a better fitting. Using the biweight location of iBeers et al.l (Il990l ). the systemic velocity of 



the M60 GC system is estimated to be Vgys = 1073^22 km s ^ (see Section [X^ . which agrees 



with the M60 recession velocit y (1056ib64 km s~^) and the biweight mean velocity (1066±45 
km s~^) for 38 GCs given by [Bridges et al.l (120061 ) within the uncertainty. We derived the 
rotation amplitudes, ^IR, 14:l~^lg km s~^ for the entire GCs, 130l^? km s"^ for the blue 



-51 



GCs, and 17ll^4g km s ^ for the red GCs. Thus, the rotation amplitude of the blue GCs is 



slightly smaller than that of the red GCs, or is consistent with that of the re d GCs within 



the un certainty. Our results based on 121 GCs are in contrast to the fact that [Bridges et al. 



(120061 ) found no rotation using only 38 GCs . To investiga t e the cause for difference of the 
rotation amplitudes between this study and [Bridges et al.l (120061 ) . we estimat e the rotation 



ampli tude using all 121 GCs by dividing G Cs into 54 GC s with in the region of [Bridges et al. 



(|2006[ ). and 67 GCs outside the region of [Bridges et al.l (|2006[ ). In the result, the rotation 
amplitude for the former sample is estimated to be 74I27 km s~^, which is much smaller 
than that for the latter sample (ilR = 221^^"^ km s^^). Therefore, it is concluded that no 



rotation found in [Bridges et al.[ (|2006[ ) is because of a small spatial coverage in their study. 



In addition, the orientation of the rotation axis (Go) is estimated to be 225°t\l for the 
entire GCs, 218°^^^ for the blue GCs, and 237°^^^ for the red GCs. The orientations of 
rotation axes for all subsamples a ppear to be simi lar, and they are closer to the photometric 
minor axis {Qphot = 15°, or 195°, [Lee et al.|[2007[ ) than the photometric major axis. Inter- 
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estingly, the position angle of NGC 4647 relative to M60 is 314° from north to east. This 
means that the line connecting M60 and NGC 4647 is nearly perpendicular to the rotation 
axes of the GC system. 

In Figure [3], we present the rotation of the GC system for the samples of different radial 
bins in order to investigate the radial variation of rotational properties. The top panel shows 
the GCs in the range of 32" < R < 533", while the lower two panels show the GCs for 
the inner region (32" < R < 200", middle panel) and the outer region (200" < R < 533", 
bottom panel). The same fitting procedure used for the results in Figure [2] was applied for 
the entire, blue, and red GCs in each radial bin. The best fit rotation curves for each sample 
are overlaid as solid (entire GCs), dashed (blue GCs), and dotted (red GCs) lines. 

It appears that the orientation of rotation axes for all GC subsamples changes slightly 
from the inner region to the outer region with a large uncertainty. The change of the rotation 
axis from the inn er region to t he outer region f or the blue GCs is seen in the cases of M87 



(ICote et al.ll200ll ) and of M49 (jCote et al.ll2003l ). The orientation of the rotation axis for the 



blue GCs in M87 drastically changes from the major axis in the inner region to the minor 
axis in the outer region. However, the blue GCs in M49 show the opposite change of rotation 
axis from the minor axis in the inner region to the major axis in the outer region. Due to the 
limited number of GCs in this study, it is difficult to make a solid conclusion on the radial 
change of the rotation axis for the M60 GC system. A larger sample of GCs is needed for 
further study on the change of the rotation axis as a function of the radius. 



3.2. Velocity Dispersion of the Globular Cluster System 

We summarize the kinematics of the M60 GC system derived in this study in Table [TJ 
Several kinematic parameters for the entire, blue, and red GCs are presented for the entire 
region (32" < R< 533"), the inner region (32" < R< 200"), and the outer region (200" < 
R < 533"). The column (1) defines the range of the projected radial distance from the 
center of M60 for each region in arcsec, and the column (2) gives the median value of the 
radial distance in arcsec. The number of GCs in each region is shown in the column (3). 
The column (4) a nd (5) represent the mean line-of-sight velocity (the biweight location of 



Beers et al.l 1X9901) a nd the velocity dispersion about this mean velocity (the biweight scale 



of iBeers et al.lll990l ). respectively. The position angle of the rotation axis and the rotation 
amplitude estimated using equation ([1]) in each region are given in the column (6) and 
(7), respectively. The column (8) gives the velocity dispersion about the best fit rotation 
curve. The column (9) gives the absolute value of the ratio of the rotation amplitude to 
the velocity dispersion about the best fit rotation curve. The uncertainties of these values 
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represent 68% (la) confidence int ervals that a re det ermined from tlie numerical bootstrap 



procedure following the method of ICote et al.l (120011 ) . The estimated velocity dispersion for 



entire 121 GCs (234^^^ km s"^) agrees with the value (256 ± 29 km s ^) of [Bridges et al.l 



(120061 ) based on 38 GCs within the uncertainty. In addition, it is found that the velocity 
dispersion about the mean velocity of the GC system for the red GCs (cTp = 258^3} km s~^) 
is margi nally larger than tha t for the blue GCs (cxp = 223l|g km s~^), which confirms the 



result of [Bridges et al.l (120061 ) with improved statistics. Interestingly, the velocity dispersion 
about the best fit rotation curve for the red GCs {<Jp,r = 240I34 km s~^) is also marginally 
larger than that for the blue GCs {ap^r = 207l}g km s~^). 

In Figure HI we plot the radial velocities of GCs with measured uncertainties against 
projected galactocentric distances. The mean radial velocities in two radial bins are overlaid 
by squares with long horizontal errorbar. The velocity dispersion about the mean velocity in 
each bin is also represented by a vertical errorbar. The mean velocities of all samples agree 
well with the systemic velocity of M60. Both of the velocity dispersions about the mean 
velocity and about the best fit rotation curve in the inner region are marginally larger than 
or comparable to those in the outer region for all three samples (see Table [1]) . 

To investigate the radial variation of velocity dispersion in detail, we present a smoothed 
radial profile of velocity dispersions about the mean radial velocity (filled symbols) and about 
the best fit rotation curve (open symbols) in Figure [51 We calculate the velocity dispersion 
of the GCs lying within a bin with fixed radial width, AR = 120" ~ 10.06 kpc as increasing 
the bin center by a fixed step width, 6R = 10" ~ 0.84 kpc. We define the radial width and 
the step width so that the number of GCs per bin exceeds 10, and the calculation stops when 
the number of GCs in a bin is less than 10. The velocity dispersions about the mean radial 
velocity of the entire and blue GCs are nearly constant for the range of radius. However, 
velocity dispersions of the red GCs are decreasing in the inner region {R < 18 kpc) and are 
increasing in the outer region {R > 18 kpc) although the variation is not large. The velocity 
dispersions about the best fit rotation curves of all three samples are not different from those 
about the mean radial velocities. 



3.3. Velocity Anisotropy of the Globular Cluster System 

If we assume the spherical symmetry of the M60 GC system, we can apply the Jeans 
equation in the absence of rotation to the dynamical analysis of the GC system. The spherical 
Jeans equation is 
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totl 



(2) 



where r is a three dimensional radial distance from the galactic center, 7T.ci(r) is a three di- 
mensional density profile of the GC system, cTrif') is a radial component of velocity dispersion, 
/3ci(r) = 1 — c'"^(r)/cj^(r) is a velocity anisotropy, G is the gravitational c onstant, and M, 



is a t otal gravitating mass contained within a sphere of radius r (e.g., iBinney fc Tremaine 
19871 ). (J0{r) is a tangential component of velocity dispersion that is equal to an azimuthal 
component of the velocity dispersion, <J(j,{r), in a spherical case. 

Several studies on the dynamics of the GC system have focused on determining the 
gravitational mass, Mtat(r), using the Jean s equation by assuraing a simple isot ropic orbit 



with /3ci(r) = fe.g.. lCohen fc Rvzhovlll997l : iMinniti et al.lll998l : IZepf et al.ll2000f ). However 



with an aid of an independent determination of the mass profile of an elliptical galaxy using 



X-ray data (e.g., iBrighenti &: Mathewd 119971: iHumphrev et al.ll200q for M60), the velocity 



aniso tropy itself can be investigated (e.g., iRomanowskv fc Kochanekll200ll : ICote et al.ll2001 



2003h . Fol lowing the analysi s of the M87 GC system bv lCote et aP J200lh and the M49 GC 
system by lCote et al.l (120031 ) . we derive first the three dimensional density profile of the GC 
system, nci(r) and the total mass profile, Mtot(?")- Comparing the velocity dispersion profile 
(VDP) calculated from the Jeans equation using those inputs with the measured VDP cTpi^R), 
we determine the velocity anisotropy of the M60 GC system. 



3.3.1. Density Profiles for the GC system 



We used the surface number density profiles of M60 GCs in iLee et al.l (120071 ). They 
derived the surface density profile of M60 GCs by combining the HST/WFPC2 archive data 
for the inner region at i? < 1.5', and the KPNO data for the outer region at i? > 1.5'. They 
determined the background levels from the mean surface number density of the point sources 
at 9' — 10' with the same range of magnitude and color as the GCs in the KPNO images: 
1.988 ± 0.363 per square arcmin for the entire GCs, 1.790 ± 0.344 per square arcmin for 
the blue GCs, and 0.199 ± 0.115 per square arcmin for the red GCs. Then they subtracted 
these background values from the original number counts to produce the radial profiles of 
the net surface number density of GCs. Since they selected GCs that are brighter than 
Ti 23.0 ma^, it is needed to correct the surface number density profile in order to account 



^This limiting maRnitu de corresp o nds t o V « 23.31 mag using the transformation relation between VI 
system and CTi system of iLee et all (j2007l ). 
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for the uncounted G Cs due to the limi ting magnitude. To calculate the correction factor, 
the equation (11) in iMcLaughlinI (119991 ) with Vum,! = oo was used on the assumption that 
the GC luminosity function of M60 has a Gaussian sh ape with a peak at ^ ~ 23.8 mag and 
a dispersion a = 1.65 mag ( Kundu fc Whitniorel 2001). It is found that the surface number 
density of the bright M60 GCs in iLee et al.l (120071 ) should be multiplied by 2.61 to derive 
the total surface number density. 

We display the total surface number density profiles, Nci{R), for the entire, blue and 
red GCs in F i gure P We fit the surface number density profile with the projection of 
Navarro et al. Jl997l . NFW) density profile, rici(r) = nn(r/b) ^(1 + r/b) ^ and with the 
projection of one of the galaxy models developed by iDehnenI (119931 ) . nci(r) = no{r / a)~'^ {1 + 
r/a)"'~'^. The surface number density profile, Nci{R), is related to the three dimensional 
density profile nci(r) as follows: 



2 / nd(r) 

'R 



r dr 



(3) 



The solid and long dashed lines represent the projected best fit curves of the NFW 
profile and of the Dehnen profile, respectively. The fitting results for the entire (E) GCs, 
blue (B) GCs, and red (R) GCs are summarized as follows: 



nfi(r) 
n^i(r) 
n^(r) 



for the NFW profile, and 



n|fi(r) 



0.61 kpc"^(r/5.96 kpc)-i(l + r/5.96 kpc)-^ 
0.32kpc"^(r/6.07kpc)-i(l +r/6.07kpc)-2 
0.46 kpc"^(r/4.99 kpc)-i(l + r/4.99 kpc)-^ 



1.42kpc"''^(r/7.97kpc)-°-29(l + r/7.97kpc)-3-^i 
0.57kpc"^(r/8.75kpc)-°-^°(l + r/8.75 kpc)^3.60 
0.53kpc~^(r/8.36kpc)-°-56(l + r/8.36 kpc)-3■^^ 



for the Dehnen profile. 



(4) 



(5) 



It is found that the scale length h of the red GCs in the NFW profile is smaller than 
that of the blue GCs, indicating that the red G Cs are more concentrated toward th e galaxy 
center than the blue GCs, as shown previously (IForbes et al.ll2004l : iLee et al.l 120071 ). 



3.3.2. Need for an Extended Dark Matter Halo in M60 



In the left panel of Figure [Tj we plotted the surface brightness p rofile of M60 derived 
from our KPNO Ti-band images (ILee et al.ll2007l ) compared to that in lPeletier et al.l (Il990l ) 
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for the i?-band photometry. We convert Ti photornetry o f iLee et aLl (120071 ) to Cousins R- 
band photometry using the relation given by lGeislerl (119961 ) . It is seen that two profiles agree 
well over the radius. 



We fit the surface brightness profile derived from the KPNO images (IKim et al.l 12006 
Lee et al. 12007 ) with the projection of three dimensional luminosity density profile used in 
Cote et al.l (120031 ). which is represented by, 



(3-7)(7-27) L 



tot 
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Ltot 


= 1.32 X ] 



2(7-4) 



(6) 



The fit yields the parameters of 7 = 0.32, Ltot = 1-32 x 10 L^,©, and a = 1.48 
kpc, and the projected best fit curve is overlaid in Figure [71 The fitted model also gives 
an effective radius of -Refr = l-'96 ~ 9.86 kpc, which is slightly lar ger than that fr om a fit 
{Res = l-'83 ~ 9.23 kpc at Ti-band) using a de Vaucouleurs law in iLee et al.l (120071 ). 



In the right panel of Figure [TJ we show a three dimensional stellar mass density profile 
Ps(r) = Toj(r), with _R-band mass-to-light ratio Tq 

Section). Thus we obtain a stellar mass profile of M60, which is represented by 



6.0 MqL^q (discussed later in this 



M,(r) 



A-Kx^ Ps{r)dx = I Atxx^ j{x)dx 



tot 



' (r/a)i/2 ■ 


2(3-7) 


■(7-27) + (r/a)^/2 


_1 + (r/a)V2_ 




1 + (r/a)V2 



(7) 



We used this stellar mass profile to determine the velocity anisotropy for the M60 stellar 
system and to test the existence of an extended dark matter halo. If we take M^o^_{r) = Ms{r) 
and substitute nci(r) by Ps(r) oc j(r), then we can compute the intrinsic radial VDP of the 
stars through the Jeans equation by assuming several i?-band mass-to-light ratios (Tq) and 
velocity anisotropies of the stellar system (/5s(r)). We therefore predict the projected VDPs 
for the stellar system from the intrinsic radial VDP calculated using equation (ITUl) . In 
Figure [HI we show the projected VDPs calculated with Tq = 6.0 MqL'^^ and /?s(r) = -1-0.6 



(radia ll y biased), which is the best fit curve for the stell ar kinematic data of [Fisher et al. 



(Il995[ ) , [de Bruyne et al.[ (l200l[ ) , and [Pinkney et al.l (120031 ) . This i?-band mass-to-light ratio 
and radially biased velocity anisotropy for the stellar syste m is similarly fou nd in the M49 



stellar system with Tq = 5.9 MqL^^ and I3s{r) = +0.3 (iCote et al.l [20031 and references 
therein). 

For the comparison, we also present the projected VDPs calculated using the same stellar 
mass profile as above, but for the GC number density profile nci(r) and for /5ci(r) = +0.99 
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(radially biased, upper long dashed lines), —99 (tangentially biased, lower long dashed lines), 
0.0 (isotropic, solid line for the NFW profile and short dashed line for the Dehnen profile). 
Interestingly, none of these models can account for the observed VDPs for the GCs at i? > 7 
kpc, indicating that mass-to-light ratio is not constant over the galactocentric distance, but 
should be increased as the distance increases. This means that there exists an extended 
dark mater halo in M60. This result is consistent with the previous findings of an extended 
dark mater halo in M6 from a radial p rofile of an increasin g mass-to-light ratio in i^-band 
jHumphrev et allbooeh and in V"-band teridges et allbooeh at 7 < i? < 22 kpc. 



3. 3. 3. X- my Mass Profile 

The total gravitating mass profile is determined using a gas temperature profile and a 
density distribution obtained from X-ray observational data on the assumption of hydrostatic 
equilibrium (neglecting magnetic pressure term): 



kT{r)r f d log p(r) ci log T(r 



tot\ 



Gfimp \ d log r 



+ 



dlogr 



where k is the Boltzmann constant, G is the gravitational constant, T(r) is a gas 
temperature at radius r, /i is the mean molecular weight (taken as 0.63 in this study), 
rrip is the proton mass, an d p(r) is a gas density at radius r. For the case of M60, 
Brighenti fc Mathewsl (1997) de rived the total mass profile using Einstein H RI observational 



1986) and ROSAT PSPC observational data of iTrinchieri et al. 



data of iTriiichieri et al 

fll997h . and iHumphrev et al.l (l2006f) derived th e total mass profile using Chandra observa- 



tional data. In addition, iRandall et al.l (120061 ) reported the density and the temperature 
distribution of M60 using XMM-Newton data. Since the derived VDP is sensitive to the 
mass profile, we estimate mass profiles using these different X-ray data set. 

In Figure [9l we display the deprojected profiles of the gas te mperature (top panel) 
and the gas number densities (midd le and botto m panels) found bv ITrinchieri et al.l (11986 



filled circles') . ITrinchieri et al.l (119971 . open circles) , iRandall et al.l (l2006l . open triangles) , and 
Humphrey et al.l (120061 . open stars) with measured errors. Having done this, we fit the tem- 
perature and gas number density data of each reference using the temperature distribution 
of T(r) = 2Tm[rm/{r + rot) + {r/rmy]~^, and the density distribution of r;,(r) = Snj(r), where 
ni{r) = nQ{i)[l + {r/ro(i)}^''*^]~"'^. T^, r^. Tot, q, nQ{i), ro(z), and p{i) are parameters of the 
fit. We fit the gas number density profile of each data for two cases of i < 1 (one component 
fit) and i < 2 (two component fit). The data in use are annotated with the associated lines. 
The results for the fit are summarized in Table [2l 
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The resulting mass profiles of M60, M(r), for several temperature and gas number 
density profiles are presented in Figure [TOl The upper panel shows the mass profiles using 
one component fit of the gas number density for various data set. The lower panel shows 
those using two c omponents of the gas num ber density for various data set and using the mass 
profile derived by [Humphrey et al.l (120061 ) in comparison. It appears that most mass profiles 
agree well. However, the mass profile d erived using the temp erature and the gas number 
density data (with one component fit) of iTrinchieri et al.l (119971 . the dotted line in the upper 
panel) deviates fr om the other profi l es in t he inner region (r < 10 kpc). In addition, the mass 
profile derived by [Humphrey et al.l (l2006l . the short dashed line in the lower panel) deviates 
from the other profiles in the region of r < 1 kpc. Comparing X-ray mass profiles with the 
stellar mass profile (heavy solid line) determined in §3.3.2[ we find that most X-ray mass 
profiles deviate significantly from the stellar mass profiles in the inner region at r < 2 kpc, 
indicating that X-ray mass profiles are not reliable at the very small radius because of angular 
resolution limit and non-equilibrium energetics. In the intermediate region at 2 < r < 10 
kpc, the X-ray mass profiles derived using two component fit of the gas number density 
(lower panel) agree to the stellar mass profiles, while those derived using one component fit 
of the gas number density (upper panel) do not. In the outer region at r > 10 kpc, all X-ray 
mass profiles are larger than the stellar mass profiles, confirming the need of an extended 
dark matter halo as shown in §3.3.21 Since the discrepancy between X-ray and stellar mass 
profiles is significant only in the inner region at r < 2 kpc, where there are no GCs, and 
X-ray mass profiles that account for the dark matter halo are similar to or larger than the 
stellar mass profile at r > 2 kpc, we conclude that the X-ray mass profiles are good enough 
to determine the velocity anisotropy of the M60 GC system for the following analysis. 



3.3.4- Determination of the Velocity Anisotropy 

We determine the velocity anisotropy of GCs as follows: (1) With the GC number 
density profile (nci(r)) of the entire, blue, and red GCs and the mass profile (Mtot(r)) in 
hand, assuming the velocity anisotropy (/5ci(r)) in prior, we derive the theoretical projected 
VDP ((Tp(i?)) and theoretical projected aperture VDP (o"ap(< R)) using the Jeans equation; 
(2) From the comparison of these calculated VDPs with measured VDPs, we determine the 
velocity anisotropy of GCs. 

We begin by deriving the theoretical projected VDPs. The equation ([2]), spherical Jeans 
equation can be solved for the radial component of velocity dispersion, ar{r): 
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Then the projected VDP, ap{R) can be derived by 

2 



ncicr^(r) 1 
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r dr 
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(10) 



where R is the projected galactocentric distance and the surface density profile, Nc\{R), is 
a projection of the three dimensional density profile nci(r). The projected aperture VDP, 
o"ap(< R), which is the velocity dispersion of all objects interior to a given projected radial 
distance R, can be computed by 
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where Rmin is the projected galactocentric distance of the innermost data point in the GC 
sample {Rmin = 2.7 kpc in this study). 

We present the measured VDP in comparison with the VDPs calculated by assuming 
several velocity anisotropies in Figures [TT] and [121 Figure [TT] shows the VDPs calculated 
using the Dehnen profile for the GC number density, while Figure [12] shows those using 
the NFW profile for the GC number density. The upper panels show the projected VDPs, 
and the lower panels show the projected aperture VDPs. The measured dispersion data 
taken from Figure [5] are shown by filled circles along with their confidence intervals. The 
projected aperture VDPs in the lower panels are plotted in the similar fashion to the case 
of the upper panel. The calculated VDPs in the left panels are obtained with the mass 
profiles derived using one component fit of the gas number density, while those in the right 
panels are obtained with the mass profiles derived using two component fit of the gas number 
density. Although it is difficult to distinguish the velocity anisotropy clearly for nearly all 
radial distances in the upper panel (the bottom panels show a more stable result), it appears 
that the entire GC system does not have an isotropic velocity ellipsoid {Pci = 0), but have a 
modest tangentially biased velocity ellipsoid [jSci < 0) for any mass profile in Figure [TT] based 
on the Dehnen profile for the GC number density. A sim ilar result can be fou nd in Figure [12] 
based on the NFW profile for the GC number density. [Bridges et al.l ( 120061 ) reported that 
the orbits of M60 GCs are close to isotropic within 100" (~8.4 kpc) and becomes tangentially 
biased beyond 100". Our results appear to be consistent with those of [Bridges et al.l (120061 ) 
as a whole, although the signature of an isotropic orb it within 100" i s wea ker in this study. 
However, we extend the radial coverage (~21 kpc) of [Bridges et al.[ (l2006l ) out to ~40 kpc 
in this study. 



In Figures [T3] and [Ml we show a similar analysis for the blue and red GCs, respectively. 
We present the results using the Dehnen profile of GC number density for the blue and red 
GCs. The results using the NFW profile are not different from those in Figures [13] and HH 



- 15 - 



Although it is not easy to draw a strong conclusion due to the small number statistics and 
complex mass profiles, we note a difference of velocity ellipsoids between the blue and red 
GCs in the projected aperture VDP (the bottom panels in Figure [TSl and UM. It appears 
that the blue GC system has a tangentially biased velocity ellipsoid with Pel < 0, while the 
red GC system has a radially biased or an isotropic velocity ellipsoid. 



4. Discussion 



4.1. Comparison with the GC Systems in Other gEs 



To date there are five giant ell iptical galaxies except for M60 that the kinematics of 



1998 
2004 



Cote et al.ll200lh. M 49 



Minniti et al. 



1998 



Zepf et al 



their GC systems was st udied: M87 (ICohen &: Rvzhoylll997l : iKissler-Patig 



200 



I 



Cote et al. 



2003h. NGC 139 9 jKissler-Patig et al 



fc Gebhardt 



19981: 



Kissler-Patig et all 19991: iRichtler et al.ll2004[) . NGC 5128 flPeng et al 



Woodlev et allboOTh . and NGC 4636 JSchuberth et al.lbood ). 



For the comparison of the kin ematic properties of GCs in those g Es, we analyze th e 



velocity data of 276 GC s in M87 (ICote et al.l l200lh. 263 GCs in M49 



435 GCs in NGC 1399 jRichtler et al.ll2004h. 210 GCs in NGC 5128 JPeng et al 



Cote et al. 



2003h . 



200^ 



and 172 GCs in NGC 4636 (ISchuberth et al.l l2006l ) using the similar method adopted in 
this study. M87, M49, and NGC 4636 as well as M60 are gEs in the Virgo cluster, and 
NGC 1399 is a gE in the Fornax cluster at the similar distance to that of the Virgo cluster. 
NGC 5128 is in the Centaurus group, being the nearest gE. Basic photometric properties of 
these galaxies are listed in Table [31 The column (1), (2), and (3) give the name of galaxy, 
the absolute magnitude in the V band, and the systemic radial velocity, respectively. The 
effective radius, the ellipticity, the position angle of photometric minor axis in degrees east 
of north, and the distance are shown in the column (4), (5), (6), and (7), respectively. The 
resulting global kinematics for the entire, blue and red GCs of each galaxy are presented 
in Table HI the mean radial velocity, velocity dispersion about the mean radial velocity, 
rotation axis, rotation amplitude, rotation-corrected velocity dispersion, and the absolute 
value of the ratio of the rotation amplitude to the velocity dispersion about the best fit sine 
curve for each GC subsample. We divide the GCs in each gE into the blue and red GCs 
using the colors that were used in the associated reference, excluding the GCs without color 
information, for the kinematic analysis. In addition, we estimate the rotation amplitude 
and the position angle of rotation axis for the GC system of NGC 5128, by fitting the mean 
velocity of GCs lying in a fixed width (60 deg) of position angle for a stable computation. 
The velocity dispersions derived in this study agree well with those derived in the associated 
reference. The velocity dispersions of GCs in M87 (cTp = 414l]^g km s^^) and NGC 1399 
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(cTp = 323]'^}3 km s ^) are found to be larger than and comparable to that of GCs in M49 
(dp = 322+}^ km s"^), respectively, although M87 and NGC 1399 are fainter than M49. This 
is due to the fact that M87 and NGC 1399 are located in the center of a galaxy cluster, while 
M49 is not. NGC 5128 shows the smallest velocity dispersion (cXp = 129~^^ km s~^) among 
our sample galaxies, implying a relatively smaller mass in spite of its high luminosity in the 
V band {My = -21.7). 

Table O lists several notable features based on the global kinematics presented in Table |H 
The strength of rotation is defined in the column (4): strong for QR/ap^r > 0.4, modest for 
0.4 > QR/ap^r > 0.2, and weak for QR/ap^r < 0.2. The rotation axis in the column (5) is 
assigned if the difference between the position angle of rotation axis and that of photometric 
major/minor axis is less than 30°. The result of velocity anisotropy for each GC system is 
taken from the literature (see column (6)). It appears that the rotation-corrected velocity 
dispersion, ap^r of the blue GCs is similar to or larger than that of the red GCs except for 
M60. This implies that the blue GC system is dynamically hotter than the red GC system 
in most gEs. However, the case of rotation (QR/ap^r) is not simple. Both of the blue and 
red GCs in M60 and M87 show a strong rotation, while those in NGC 1399 show a weak 
rotation. The red GCs in NGC 4636 and NGC 5128 show a slightly stronger rotation than 
the blue GCs. However, the red GCs in M49 show a weak rotation, while the blue GCs 
show a modest rotation, wh ich is consistent with the prediction of the merger formation 



model (j Ashman fc Zepflll992l ). In addition, if we consider the position angle of rotation axis 
together, the story becomes more complicated. In velocity anisotropy, it appears that the 
entire sample of the GC system in three gEs (M87, M49, and NGC 1399) has an isotropic 
velocity ellipsoid and two subsamples have different velocity ellipsoids in two gEs (M60 and 
M87). It is needed to determine the velocity anisotropy of the GC systems in more gEs 
including NGC 5128 and NGC 4636. 

To examine the general kinematic properties of the GC systems in gEs, we plot the 
rotation-corrected velocity dispersions in gEs against the projected galactocentric distances 
in Figure [151 The rotation-corrected velocity dispersion is normalized by that of the entire 
GCs in each gE. The projected galactocentric distance is normalized by the effective radius 
of each gE. The entire and blue GCs do not show significant change of the velocity dispersion 
over the whole region of a galaxy. However, the velocity dispersion of the red GCs in the 
inner region {R < 2Rcs) is, in the mean, marginally larger than that in the outer region 
{R > 2Rcs)- This implies that red GC system in the inner region may be dynamically hotter 
than that in the outer region. In Figure dni the absolute value of the ratio of the rotation 
amplitude to the velocity dispersion is plotted as a function of the projected galactocentric 
distance. For the red GCs, the ratio of the rotation amplitude to the velocity dispersion in 
the outer region appears to be marginally larger than that in the inner region if we neglect 
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the point of NGC 4636 (filled pentagon), while it does not change with the distance for the 
entire and blue GCs. 

In summary, differences in kinematic properties among GC subsamples appear to exist, 
although it depends on the galaxy. The blue GC system appears to be dynamically similar 
to or hotter than the red GC system, while the rotation of the GC system is not negligible. 
The entire sample of each GC system appears to have an isotropic velocity ellipsoid, while 
the subsamples do not have uniform velocity ellipsoids. The kinematic properties of M60 and 
M87 are similar except for the velocity dispersion of subsamples and the velocity anisotropy 
of the entire GC system, while other galaxies have diverse kinematic properties. For the red 
GCs, the velocity dispersion of the inner region is marginally larger than that of the outer 
region, while the rotation of the outer region appears to be more significant than that of the 
inner region. 



4.2. Formation Models of Globular Clusters 



The kinematic properties of the GC system in M60 compared to those in other gEs are 
useful to test the kinematic predictions of formation models of the GCs in gE. A summary of 
sever a l model d e scriptions and predictions can be found in several literature (jRhode fc Zepi 
200ll : Q I2OO3I : iRichtler et all I2OO4I : IWest et al.l bood : ISrodie fc Stradeil I2OO6I). Classical 
formation models can be divided into four broad categories: the monolithic collapse model, 
the major merger model, the multiphase dissipational collapse model, and the dissipationless 
accretion model. 

The monolithic collapse model describes that an elliptical galaxy and its GCs are formed 
through the collapse of an isolated massive ga s cloud or protogalaxy at high redshift (ILarson 
I975I : ICarlberg||l984l : lArimoto fc Yoshiil 119871 ). In this model, the color distribution of GCs 
shows a smooth shape with a single peak , and the rotation of GCs can be generated by tidal 



torques from companions (jPeeblesI Il969l ). Although t his model can explain some observa- 
tional properties of elliptical galaxies successfully (see IChiosi fc Carraroll2002h. the bimodal 



color distribution of the GC system in many gEs (e.g., Leel 20031 : Peng et al. 2006) make 



it hard to accept this model. In addition, the strong rotation of the GC system as seen 
in M60 and M87 is not expected from the collapse of a single, isolated protogalactic cloud. 
Moreover, several GC systems in gEs show a globally isotropic velocity ellipsoid that needs 
some kinds of relaxation processes that the monolithic collapse model can not account for. 



The major merger mo del suggests that elliptical galaxies are formed by a merger of 
two or more disk galaxies (IToomrd 119771 : lAshman fc Zepj 1 19921 : IZepf et al.l 12000) • In this 
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model, younger, spatially concentrated, red GCs are formed during the merger, while spa - 



tially extended, blue GCs come from the halos of the disk galaxies (e.g., iBekki et al.ll2002l ). 
As a result, the color distribution of GCs is expected to be bimodal. This model predicts 
that the newly formed red GCs show little rotation compared to the blue GCs since the 
angular momentum would be transported to the outer region during the merging process. 
This model has received particular attention since it could explain several photo metric prop 



erties of the GC system in gEs and little rotation of the red GC system in M49 (jZepf et al. 



200(1 1. However, contrary to the case of M49, the red GC system in M60, M87, NGC 5128, 
and NGC 4636 show a strong or modest ro t ation . From the simulation of dissipationless 
major mergers of spiral galaxies, iBekki et al.l (120051 ) found that both pre-existing metal-poor 
clusters (MFCs) and metal-rich clusters (MRCs) obtain significant amounts of rotation be- 
yond ~ 10 kpc, regardless of the orbital configuration of the merging galaxies. However, 
both the blue and red GC systems i n NG C 1399 show weak rotation, which is not consis- 
tent with the result of iBekki et al.l (120051 ). Therefore, it is needed to explain what makes 
the complex rotational properties of the GC systems in gEs. In the dissipationless major 
merger model, most GCs are formed before a last major merger event. There are several 
observational results based on the spectroscopy of the GCs in gEs consistent with this: the 



blue GCs an d red GCs show a small a ge difference and ar e both old : M60 (IFierce et al. 



2006h . M87 (ICohen et a: 



1399 (IKissler-Fatig et al. 



1998; 



1998 ). M49 (IBeaslevet al.l l2000l: ICohen et all 120031 ). and NGC 



Forbes et al.ll200ll ). Then, what we observe today is the GC 



system that was affected by the orbital mixing caused by the last major merger event, and 
the observed kinematics can not reflect the circumstance when the GC system was formed. 
Therefore, it is noted tha t it is difficult to trace the orbital history of the GC system from 
the observed kinematics (IKissler-Fatig fc Gebhardtl 119981 ). and to test the validity of the 
dissipationless merger scenario. 



Forbes et al.l (119971 ) proposed the multiphase dissipational collapse model that ellipticals 



form their GCs in distinct star formation phases through a dissipational collapse. In addition, 
there is a capture of additional GCs by tidal effects from neighboring galaxies or the accretion 
of dwarf galaxies. Since the blue GCs are formed in the first star formation phase and the 
red GCs are formed in the subsequent star formation phase after the gas in the galaxy is 
self-enriched, the color distribution of GCs is expected to be bimodal. This model predicts 
that the blue GC system shows no rotation and a high velocity dispersion, while the red GC 
system shows some rotation depending on the degree of dissipation in the collapse. As the 
galactocentric distance incre ases up to 90", the stell ar velocity dispersion in M60 approaches 
to the value of 200 km s~^ (jde Bruyne et al.ll200ll ) that is comparable to or smaller than 



the ve locity dispersion of the blue GCs. This might support this model, as iForbes et al. 



( 119971 ) pointed out in the case of NGC 1399 and M87. However, this model disagrees with 
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the results that blue GC systems in several gEs including M60 in Table [5] show larger than 
or comparable rotation amplitude to that of the red GC system. 



Cote et al.l (Il998l ) proposed the dissipationless accretion model; the red GCs formed in 
a dissipational collapse, while the blue GCs were subsequently captured from other galaxies 
through mergers or tidal stripping. This model predicts a bimodal color distribution of GCs 
in gEs and a uniform color distribution in dwarf ellipticals. However, there is an example 
of low- luminosity elliptical galaxy (e.g., NGC 1427) that shows a bimodal color distribution 
( iForte et al.ll200ll ). Since the blue GCs are captured f r om ot her galaxies, they are expected 
to show ex t ended spatial distribution. iRichtler et al.l (120041 ) argued that if the scenario of 
Cote et al.l (Il998l ) is correct, then the blue GCs are expected to have radially biased orbits 
rather than the isotropic or tangentially biased orbits, and are also expected to show no 
rotation. Although the kinematic data for the outer GCs are not complete to date, current 
data for the blue GCs in Table [5] show negligible rotations and no signs of radially biased 
orbits, which is not consistent with the scenario of ICote et al.l (119981 ). 



Since the classical models do not give quantitative predictions concerning the kinematic 
properties of the GC system, it is instructive to compare the observational results with those 
in the recent numerical simulations. Although there are sonie num erical simulations of the 
GC system focusin g on the color distribut ion (IBeasley et al.ll2002l ). the spatial distribution 
(iMoore et al.l 120 06: Be kki fc Forbes! l2006l ). and the mass-metallicity relation for the blue 
GCs (IBekki et al.ii2007l ). there are few simu lation results that can be compared directly w ith 
observational results summarized in Table (IBekki et al.ll2005l : iKravtsov fc Gnedinll2005l ) 



In a pioneering work, iBekki et al.l (120051 ) numerically investigated the kinematics of the 
GC system in E/SO galaxies formed from a dissipationless merging of spiral galaxies. They 
presented the kinematic properties such as rotation and velocity dispersion of the pre-existing 
MFCs and MRCs for several merger configurations (e.g., pair and multiple mergers). For 
the rotation, both MFCs and MRCs show larger rotation amplitudes in the outer region 
{R > 2Res) than those in inner region [R ~ -Reff), regardless of the merger configuration. 
Interestingly, observational data in Figure [16] show that the rotation amplitudes for the blue 
GCs do not change as the galactocentric distance increases. However, those for the red GCs 
increase marginally from the inner region to the outer region if we neglect the point of NGC 
4636 (filled pentagon), as expected in the simulation. For the velocity dispersion, the MFCs 
show slightly larger the central velocity dispersion than the MRCs, indicating that the MFCs 
are dynamically hotter than the MRCs. Moreover, the VDFs for both MFCs and MRCs are 
decreasing as the galactocentric distance increases in all major merger models, while those 
are sometimes very flat in multiple merger models. The kinematic data in Figure [15] show 
that the VDFs might be different depending on the galaxy (e.g., flat for M49 and increasing 
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for M87), and the difference of VDPs between blue and red GCs appears to exist. Although 
the latter is not expected in the simulation, the former is consistent with the simulation, 
implying different merger histories depending on the galaxies. 



Kravtsov fc Gnediru (120051 ) studied the formation of the GCs in a Milky Way-size galaxy 
using a gas dynamics cosmological simulation with an adaptive refinement tree code. How- 



ever, t hey had to stop t 



L ater, 
in 



l e simu lation at 2; ~ 3.3 due to the limited computational resources. 



Gnedin fc Prieto 



(120061 ). using a separate coUisionless A^-body simulation described 
Kravtsov et al.l (j2004j ) up to 2; = 0, calculated the GC orbits for the similar galactic sys- 



tem. They found that, at present, the GC orbits are isotropic in the inner 50 kpc region 
from the galactic center, but radial in the outer region. Although it is difficult to exam the 
observational data at large distance (> 50 kpc) because of insufficient GC samples at that 
region, the isotropic orbit in the inner region is roughly consistent with the kinematic data 
of gEs (see Table [5]) . 



On the other hand, IVesperini et al.l (120031 ) modeled a dynamical evolution of the M87 
GC system with different initial conditions for the mass function, the spatial distribution, 
and the v elocity ani s otropy of th e GC system, following the evolution of a model GC system 
given by IVesperinil (120001 . I2OOII ). In their result, for a two-slope power-law initial mass 
function of the GCs corresponding to the mass function of the old clusters, a flat radial 
profile for the mean mass of the GCs obtained from the observation can be acquired from 
the simulation with any velocity anisotropy of the GC system. However, for a power-law 
initial mass function of the GCs that is regarded as that of the young clusters, the observed 
fiat radial profile for the mean mass of the GCs can be obtained from the simulation only with 
a strong radial anisotropy of the GC system. This initial radial anisotropy is much more 
radially biased than the observed one in the GC systems of gEs (having mostly isotropic 
velocity ellipsoids as seen in Tabled]), raising a problem for a scenario that the young GC 
system in mergers like the Antennae may evolve into the GC system seen in gEs. 

By viewing the formation of the M60 GC system in the context of kinematic predictions 
of above models, the strong r otatio ns of the blue and red GC system are consistent with 
the simulation of iBekki et al.l (120051 ). However, the larger velocity dispersion of the red GCs 
compared to the blue GCs in M60 is not consistent with the predictions of several forma- 
tion models. M60 might be a particular case among gEs due to a possible interaction with 
the c ompanion galaxy, NGC 4647. While they were regarded as a non-interacting system 



[e.g. 



Sandage fc Bedkelll994l ). recent observational results suggest evidences of current in- 



teraction between two : (1) the morphological structure of NGC 4647 is clearly asymmetric 



( iKoopman et al.l 1200 ll ): (2) the stellar kinematic study shows that the inner region of M60 
has a strong rotational support compared to other gEs, and has an asymmetric rotation 
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curve (jPinkney et al.ll2003l : Ide Bruyne et al. 



2001 



north-eastern edge of M60 is seen (IRandall et a. 



associations in NGC 4647 are found (iLee et al 



; (3) an X-ray filament that extends to the 



20061) ■ (4) youn g luminous star clus ters or 



20071). However, IPierce et al.l (120061 ) found 



no obvious signs of a recent starburst, interaction or merger by estimating the ages and 
metallicities from the spectra of M60 GCs. Thus, it appears that the interaction between 
two started very recently, and di d not affect the old M60 GC system significantly. The spa- 
tial coverage of observed GCs in iPierce et al.l (120061 ) is not enough to study the interacting 
region between M60 and NGC 4647 (see Figured]). Therefore, to understand the formation 
of M60 GC system in terms of an interaction between galaxies, it is important to obtain 
a large spectroscopic sample of GCs at larger galactocentric distance with high S/N ratio 
enough to determine the age and metallicity. 



5. Summary 

Using the photometric and spectroscopic database of 121 GCs (83 blue GCs and 38 red 
GCs) in the gE M60 (NGC 4649) in the Virgo cluster, we have investigated the kinematics 
of the GC system of this galaxy. Our primary results are summarized below: 



1. Similar to the case of M87 GCs (jCote et al.ll200ll ). the entire, blue and red GC sub- 



samples of M60 show significant overall rotations. The rotation axes are nearly per- 
pendicular to the line connecting M60 and its companion NGC 4647. 

Both of the velocity dispersion about the mean velocity and about the best fit rotation 
curve of the red GCs are marginally larger than those of the blue GCs. This implies 
that the red GC system might be dynamically hotter than the blue GC system unlike 
the GC systems in other gEs. 

Comparison of observed stellar and GC velocity dispersion profiles with those calcu- 
lated from the stellar mass profile showed that the mass-to-light ratio is not constant, 
but should be increased as the galactocentric distance increases, indicating the exis- 
tence of an extended dark matter halo in M60. 

Using the X-ray mass profile, the number density distribution of GCs, and the observed 
VDP of GCs, we have determined the velocity ellipsoids of the M60 GC system. The 
entire GC system in M60 appears to have a tangentially biased velocity ellipsoid unlike 
the GC systems in other gEs. Two subsamples have a different velocity anisotropy: 
the blue GCs show a modest tangentially biased velocity ellipsoid, while the red GCs 
show a modest radially biased or an isotropic velocity ellipsoid. 
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5. We have compared the kinematics of the M60 GC system in this study with the results 
of other GC systems in gEs. As a whole, the rotation-corrected velocity dispersion 
of the blue GCs is similar to or larger than that of the red GCs unlike the M60 
GC system. The rotation of the GC system in gEs is not negligible in contrast to the 
traditional view, though the details of rotation amplitudes and rotation axes need to be 
investigated further. The entire sample of each GC system appears to have an isotropic 
velocity ellipsoid, while the subsamples do not show unified velocity anisotropy. 

In conclusion, the GC systems in gEs have common kinematic properties such as velocity 
dispersion, while the rotation and the velocity anisotropy show diverse results. These kine- 
matic properties are not fully explained by any current models. We need extensive kinematic 
studies of other GC systems over larger galactocentric distances with various environmental 

effects to understand the diversity of the kinematics, and the formation and evolution of the 
GC system in gE. Moreover, it is desirable to have more elaborate model predictions for the 
kinematics of the GC system in gE including the velocity anisotropy. 
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Fig. 1. — Spatial distribution of M60 GCs with measure d velocities. The bl ue GCs measured 
in Paper I are represented by circles and those from iPierce et al.l (120061 ) are represented 
by pentagon s . Th e squares and triangles indicate the red GCs measured in Paper I and 
Pierce et al.l (120061 ) . respectively. The GCs with larger velocities than the systemic velocity 
(vgai = 1056 km s~^) of M60 are plotted by open symbols, while those with smaller velocities 
than the systemic velocity of M60 by filled symbols. The symbol size is proportional to the 
velocity deviation. M60 and its companion galaxy, NGC 4647, are shown by large dotted 
D25 ellipses. The photometric major and minor axes of M60 are represented by the dashed 
lines. 
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Fig. 2. — Radial velocities versus position angles for entire 121 GCs (top panel), 83 blue 
GCs (middle panel), and 38 red GCs (bottom panel). The solid curve represents the best 
fit rotation curve from Table [Tj and the dot-dashed horizontal line indicates the systemic 
velocity of M60. The photometric minor axis of M60 is represented by the vertical arrows 
(e = 15° and 195°). 
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Fig. 3. — Radial velocities versus position angles for the GCs in the range of 32" < R < 533" 
(top panel), 32" < R< 200" (middle panel), and 200" < R< 533" (bottom panel). Filled 
circles indicate the blue GCs, while open circles the red GCs. The best fit rotation curves 
for the entire (solid curves), blue (dashed curves) and red (dotted curves) GCs within each 
region are overlaid. The dot-dashed horizontal line indicates the systemic velocity of M60, 
and the vertical arrows mark the position angle of photometric minor axis of M60. 
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Fig. 4. — Radial velocities versus projected galactocentric distances for the entire (top panel), 
blue (middle panel), and red (bottom panel) GCs. Large open squares indicate the mean 
radial velocities of GCs in the radial bins that are represented by long horizontal errorbars. 
Their vertical errorbars denote the velocity dispersions of GCs in the radial bins. The dot- 
dashed horizontal line indicates the systemic velocity of M60. 
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Fig. 5. — Radial profiles of velocity dispersion for the entire (top panel), blue (middle panel), 
and red (bottom panel) GCs. Filled circles indicate the velocity dispersion about the mean 
GC velocity (dp) at each point, while open circles the velocity dispersion about the best fit 
rotation curve (cTp,.) at the same point. The dispersion is calculated using the GCs within 
the moving radial bin (width of 2' ~ 9.77 kpc) that is represented by a horizontal errorbar 
in top panel. The dotted and dot-dashed lines denote 68% and 95% confidence intervals on 
the calculation of velocity dispersion, respectively. The dashed horizontal fine indicates the 
global value of velocity dispersion of GCs in each panel. In the bottom panel, large filled 
and open squares, respectively, represent the velocity dispersion about the mean GC velocity 
((7p) and that about the best fit rotation curve {(Jp^r) for the red GCs beyond i? ~ 25.1 kpc. 
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Fig. 6. — Projected number density profiles for tlie entire (top panel), blue (middle panel), 
and red (bottom panel) GC candidates. Filled circles represent the GC candidates from 
HST/WFPC2 a rchive, while open squares the GC candidates from the KPNO CTi images 
( iLee et al.l 120071 ). The solid line and the dashed line in each panel indicate the projected 
best fits using the NFW density profile and the Dehnen density profile, respectively, for each 
sample. 
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Fig. 7. — Left: f?-ban d surface photometry of M60 deri ved from KPNO ima ges (iLee et al. 
20071 : iKim et al.l |2006| . open circles) compared to that in iPeletier et al.l (jl990l . filled circles) . 
Dashed line indicates a projected best fit using eq. ([6]). Right Three dimensional stellar mass 
density profile using the best fit model in the left panel with a constant -R-band mass-to-light 
ratio of To = 6.0 MqL^q. 
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Fig. 8. — yPPs for the stars and the GCs. Stellar V DPs are from iFisher et al.l (119951 . filled 
triangles) . Ide Bruyne et al.l (1200 ll . open squares), and lPinkney et al.l (l2003l . open stars), and 
the GC VDPs shown by filled and open circles with associated, dotted and dot-dashed hues 
are from Fig. O Dotted curve represents the stellar VDP calculated using the stellar mass 
model in Fig. [7] with a constant stellar mass-to-light ratio of Tq = 6.0 MqL^q and a 
stellar velocity anisotropy of /3s = 0.3. Other lines indicate the VDPs calculated using the 
same stellar mass model as above with several GC density profiles and velocity anisotropies: 
Dehnen and /3d = 0.0 (short dashed line), NFW and (3^ = 0.0 (solid line), NFW and 
/3ci = +0.99, -99 (long dashed lines). 
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Fig. 9. — Radial profiles of the gas temperature (a), the gas number density with one 
component fit (b), and the gas number density with two component fit (c ) . Th e dif- 
ferent s ymbols represent the d ata from differe nt literature: iTrinchieri et al.l ( jl997l . open 
circles), Trinchieri et al. ( 1986 . filled circles), Randall et al. (j2006 . open triangles), and 
Humphrey et al. open stars). The dotted, dot-dashed and solid lines in each panel 



indicate the best fit curves using the annotated data, while th e long dashed lines in (a) and 
(c) denote the profiles derived by lBrighenti &: Mathews! (119971 ) . 
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Fig. 10. — Total mass profiles of M60 derived (a) using the gas number density profile with 
one component fit, and (b) using the gas number density profile with two component fit. 
The solid, dot-dashed and dotted lines indicate the mass profiles derived in this study using 
the annotated data, while the l ong dashed line and the sho rt d ashed line in (b) denote 
the M60 mass profiles derived by lBrighenti &: Mathewd (119971 ) and [Humphrey et al.l (120061 ) . 
respectively. A heavy solid line represents the stellar mass profile using a constant i?-band 
mass-to-light ratio of Tq = 6.0 MqL~^^q derived in §3.3.21 
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Fig. 11. — VDPs (a, c) and the aperture VDPs (b, d) for the entire GCs. Filled circles in (a, 
c) represent the measured VDP shown in Fig. O and those in (b, d) denote the measured 
aperture VDP. Associated, dotted and dot-dashed lines represent 68% and 95% confidence 
intervals on the calculation of velocity dispersion, respectively. Three smoothly curved lines, 
from a radially biased velocity anisotropy to the tangentially biased velocity anisotropy (from 
top to bottom, Pci= 0.99, 0, and —99), represent the VDPs calculated using the GC density 
profile of Dehnen with one component fit of the gas number density profile (a, b) and with 
two component fit of the gas number density profile (c, d). 
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Fig. 12. — Same as Figure [TTl but for the VDPs calculated using the GC density profile of 
NFW. 
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Fig. 13. — Same as Figure [TTl but for the blue GCs. 
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Fig. 14. — Same as Figure [TT], but for the red GCs. In the top panels, large filled square 
represents the measured velocity dispersion about the mean GC velocity (dp) beyond R ~ 
25.1 kpc. 
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Fig. 15. — Rotation-corrected velocity dispersions in gEs versus the projected galactocentric 
distances for the entire (top panel), blue (middle panel), and red (bottom panel) GCs. The 
rotation-corrected velocity dispersion is normalized with respect to that of the entire GCs 
in each gE. The projected galactocentric distance is normalized with respect to the effective 
radius of each gE. Open symbols indicate the dispersions in the inner region of each gE, 
while filled symbols those in the outer region. 
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Fig. 16. — Absolute values of the ratio of rotation amplitude to velocity dispersion versus 
the projected galactocentric distances for the entire (top panel), blue (middle panel), and 
red (bottom panel) GCs. Open symbols indicate the ratios in the inner region of each gE, 
while filled symbols those in the outer region. 
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Table 1. Kinematics of the M60 Globular Cluster System 



R 


{R) 


N Vp ap Go 




nR/ Cp^y 


(arcsec) 


(arcsec) 


(km s~^) (km s~^) (deg) 


(km s~^) (km s~^) 





Entire GCs: 121 Clusters with 1.0 < (C - Ti) < 2.4 



32- 


-533 


218 


121 


10731^^ 


2341^1 




1411^° 


217li^ 


0.651°!^ 


32- 


-200 


116 


60 


10461^1 




212^25 


1 Q7+83 


2301^^ 


0.60l°j° 


200- 


-533 


319 


61 


10981^^ 




24lli^ 


1561^^ 




0.78l°J° 



Blue GCs: 83 Clusters with 1.0 < (C - Ti) < 1.7 



32- 


-533 


228 


83 


1086li 


223+}^ 


2181^^ 


130+^? 


2071}^ 


fiS+°-^^ 

U-O<J-0.30 


32- 


-200 


121 


43 


10761^1 


2301^^ 


2161^ 


13911^ 


1991^1 


7n+o-67 

•J- ' U_o.45 


290- 


-533 


343 


40 


10971^^ 


2161^^ 


991 "'"^^ 
^^-•--26 


1211^^ 


2121^? 


U.d ( _o.24 



Red GCs: 38 Clusters with 1.7 < (C - Ti) < 2.4 



32- 


-533 


197 


38 


10401^^ 


2581^} 


237l}« 


171+58 


2401^2 


71+0-30 


32- 


-200 


102 


17 


9611^1 


282^25 


1931^? 


1391^4 


2811^^ 


n 40+0-38 
'-'•^^-0.24 


200- 


-533 


274 


21 


iiool^^ 


2191^1 


2491^1 


2181^1 


18611^ 


1 17+0.58 
-0.62 
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Table 2. Gas Temperature and Number Density 



Data 




T(r) 




n 


(r) (i < 1) 






n(r) (i < 2) 










Tot 


q 


noil) 


r-o(l) 


P(l) 


no(l), no(2) 


ro(l), ro(2) 


P{1), P{2) 




(10'^ K) 


(kpc) 


(kpc) 




(cm~3) 


(kpc) 




(cm~"^) 


(kpc) 




T97 ° 


1.16 


37.32 


23.03 


-0.003 


0.02 


2.60 


1.45 








T97+T86 ' 










0.09 


0.70 


1.41 


0.07, 0.0009 


0.97, 16.83 


1.69, 2.22 


R06 " 


7.47 


19.49 


17.34 


0.19 


0.14 


0.81 


1.65 


0.06, 0.05 


0.81, 1.61 


3.08, 1.71 


BM97 


0.9 


24.28 


24.28 


0.0 








0.1, 0.0014 


0.9, 18.21 


1.8, 3.0 



iTrinchieri et al.l l|l997l V 

iTrinchieri et al.l l|l997h and iTrinchieri et al.l (|l986h . 
'^ iRandall et al.l l|2006l '). 

Values taken from ISrighenti &: Mathew 



Table 3. Giant Elliptical Galaxy Samples 



Galaxy 


My " ^sys ^ RcS 


, d PA . c 
c J. Jiminor 


Distance ^ 




(km s^"*^) (kpc) 


(deg) 


(Mpc) 



M60 


-22.44 


1056 


9.23 


0.21 


15 


17.3 


M87 


-22.62 


1307 


7.66 


0.125 


69 


17.2 


M49 


-22.83 


997 


9.97 


0.175 


65 


17.1 


NGC 1399 


-21.95 


1442 


14.55 


0.099 


20 


20.0 


NGC 5128 


-21.66 


541 


6.02 


0.224 


125 


4.2 


NGC 4636 


-21.43 


906 


6.37 


0.256 


58 


14.7 



Absolu te V magnitude base d on Bt, {B — V)t (jde Vaucouleurs et al. 



199ll ). Ay (jSchlegel et al.lll998l ). and distance adopted in this study. 



M60 (This study), M87, M49 fISmith et al 



(Richtler et al. 




004 




fSchuberth et al 


2006) 



2000), NGC 1399 



Position angle o f the minor axis. M 60 flLee et all 120071) . NGC 5128 



( Dufour et al. 1979), and other galaxies ( Kim et al. 20061 ). 



^ M60 dLee et al.ll2007h NG C 5128 (Ide Vaucouleurs et all Il99lh and 
other galaxies (IKim et al.l 120061 ). 



^ M60, M87, and M49 flMei et al.l 120071 ). and other galaxies 
jTonrv et allboOlh . 
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Table 4. Global Kinematic Properties of GCs in gEs 



Galaxy 


GC 


N 


Vp 

(km s~^) 


(Tp 

(km s~^) 


©0 

(deg) 


(km s~ ) 


(km s~^) 




M60 


EGG 
BGG 
RGG 


121 

83 
38 


1073ti 
1040+!^ 


2341}! 
2231}^ 
258+1} 


2251}! 
2I8I23 
237+}^ 


14111^ 
1301^? 

171+4fi 

— 4D 


2171}^ 
2071}^ 
240+15 


0.65l°:i 
0.63l°J^ 
0.71+°-^° 

— U.zy 


M87 


EGG 
BGC 
RGG 


276 
158 
118 


133311 
134111 
1324+0? 


4141}^ 
4251^^ 


681}^ 
591}^ 
79+}^ 

' ^—17 


1811L' 
165+^o 


3991}^ 
414+22 

3801^^ 


44+0-16 

'-'•^^-0.13 

43+° }^ 
" "— 0.12 


M49 


EGG 
BGG 


263 

159 
1 n/i 


9731?^ 
954lp 
yyy_25 


3221}^ 
3521^^ 


1021^6 

1 so+53 

■•■0^-50 


54li 

921^^ 
-1-1+79 

^^-83 


3211}^ 

3491^^ 
970+19 


17+0-17 

rv 97+0.22 
'-'■^'-0.12 

n 04+0-29 
U-U4_o.3o 


NGC 1399 


EGG 
BGC 
RGG 


435 
216 
219 


14451^6 
14391}? 


323l}J 

359 

2851}^ 


3071^° 
2611^^ 

Q+40 
'-'-36 


311^^ 
691^^ 

461^^ 


3261}^ 
3641^^ 
2881}^ 


'-'•-'-'-'-0.15 

19+0-20 

'-'•-'-^-0.09 

16+°-^^ 

'-'•-'-"-'-0.15 


NGC 5128 


EGC 
BGC 
RGG 


210 
127 
83 


5361^ 

526+}} 

5521}^ 


1291^ 
1261^ 
133l?i 


184li 
1681^? 
1911^° 


301}^ 

oc:+22 
^^^-35 

47lg 


1291^ 
1291^ 
1321?! 


'-'•^'-'-0.13 

19+0-18 

'-'•-'-^-0.29 

0-36l°i^ 


NGC 4636 


EGC 
BGC 
RGG 


172 

96 
76 


8991}^ 
90lli 
896l^J 


2071}? 

2071}^ 
2081}^ 


2751^^ 

234l| 
2891^° 


90+35 

-^-'-53 

501^? 


207+}° 
2071}^ 
208l}2 


14+0-18 

'-'•-'-^-0.08 

09+°-^^ 
24+°-23 

U.Z4_Q^2 
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Table 5. Summary of Global Kinematic Properties of GCs in gEs 



Galaxy 


Velocity Dispersion 


GC 


Rotation 


Rotation Axis 


Pel 


M60 


(jp^r{BGC) <ap,r{RGC) 


EGC 
BGC 
RGC 


Strong 
Strong 
Strong 


Minor Axis 
Minor Axis 
None 


modest tangential 
modest tangential 
modest radial/isotropic 


M87 


ap,r{BGC) > ap,r{RGC) 


EGC 
BGC 
RGC 


Strong 
Strong 
Strong 


Minor Axis 
Minor Axis 
Minor Axis 


isotropic 
modest tangential 
modest radial 


M49 


ap,riBGC) > Up^riRGC) 


EGC 
BGC 
RGC 


Weak 
Modest 
Weak 


None 
None 
Major Axis 


isotropic 
closely isotropic 
closely isotropic 


NGC 1399 


ap,r{BGC) > ap,r{RGC) 


EGC 
BGC 
RGC 


Weak 
Weak 
Weak 


Major Axis 
Major Axis 
Minor Axis 


closely isotropic 
closely isotropic 
closely isotropic 


NGC 5128 


ap^r{BGC) ~ ap,r{RGC) 


EGC 
BGC 
RGC 


Modest 
Weak 
Modest 


None 
None 
Major Axis 




NGC 4636 


ap{BGC) ~ ap{RGC) 


EGC 
BGC 
RGC 


Weak 
Weak 
Modest 


None 
Minor Axis 
None 


modest tangential 



M60 (This s tudy), M87 (ICote et al.ll2nnih . M49 (ICote et al, 
and NGC 4636 (jSchuberth et al.ll20n6h . 



20031 ). NGC 1399 (|Richtler et all 120041 ). 



